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Abstract: We present a simple super symmetric axion model that can naturally explain 
the origin of the Higgs fj, and Bfi terms in gauge mediation while solving the strong CP 
problem. To stabilize the Peccei-Quinn scale, we consider mixing between the messenger 
fields that communicate the supersymmetry and Peccei-Quinn symmetry breaking to the 
visible sector. Such mixing leads to the radiative stabilization of the Peccei-Quinn scale. In 
the model, a Higgs coupling to the axion superfield generates the B parameter at the soft 
mass scale while a small ^ term is induced after the Peccei-Quinn symmetry breaking. We 
also explore the phenomenological and cosmological aspects of the model, which crucially 
depend on the saxion and axino interactions with the ordinary particles induced by the 
Higgs coupling to the axion superfield. 
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1. Introduction 

Gauge mediation of supersymmetry (SUSY) breaking |jl], ^ is an attractive mechanism to 
generate soft terms in the minimal supersymmetric Standard Model (MSSM). In particular, 
since the transmission occurs through the gauge interactions, models of gauge mediation 
solve the SUSY flavor and CP problems. However, gauge mediation has difficulty in ex- 
plaining the origin of the Higgs // and terms, and thus requires an extension. Since 
the fi term breaks the Peccei Quinn (PQ) symmetry Q, an interesting possibility is that 
the presence of a /i term has the same origin as the invisible axion ||5| solving the strong 
CP problem ^, |9[. The size of the Bfi term is then determined by how the saxion, 

the scalar partner of the axion, is stabilized. A potential for the saxion is generated only 
after SUSY breaking because the PQ symmetry makes the scalar potential flat along the 
saxion direction in the supersymmetric limit. This indicates that it is non-trivial to stabi- 
lize the saxion, which is a gauge singlet, within gauge mediation where soft terms receive 
contribution proportional to the gauge couplings. 
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In this paper, we consider a simple axion model in the framework of gauge mediation 
that provides a natural solution both to the jJi/Bji problem and the strong CP problem. The 
model contains matter fields that communicate SUSY breaking from a hidden sector to the 
MSSM, and also those that transmit the PQ symmetry breaking. They are charged under 
the Standard Model (SM) gauge groups, and become massive either by directly coupling 
to the hidden sector SUSY breaking field or by coupling to the axion superfield. The 
most important property of the model is that there is mixing between these two classes of 
messengers. It is through the mixing that the saxion feels SUSY breaking at the loop level 
and acquires soft mass comparable to those of the gauge-charged sparticles. As a result, 
the saxion is radiatively stabilized at a scale below or around the scale of gauge mediation. 
In this model, the /v, term arises from an appropriate coupling of the Higgs doublets to the 
axion supermultiplet either in the superpotential or in the Kahler potential. Remarkably, 
the term is then generated at the correct mass scale, thanks to the SUSY breaking in 
the axion supermultiplet induced by the mixing between the two messenger sectors. 

The MSSM soft terms receive negligible threshold corrections at the PQ scale, but their 
renormalization group (RG) evolutions are affected by the PQ messengers. In particular, 
if the saxion has a vacuum expectation value rather close to the scale of gauge mediation, 
MSSM scalar masses can receive a sizable contribution through the hypercharge trace 
term because the PQ messenger scalars acquire additional soft masses due to the mixing. 
This contribution can make the stau lighter than the gauginos. In the model, the lightest 
superparticle (LSP) is given either by the axino, the fermionic partner of the axion, or by the 
gravitino depending on the scale of gauge mediation. The ordinary sparticles dominantly 
decay into axinos, not into gravitinos, through the interactions suppressed by the PQ scale. 
Meanwhile, the saxion properties are constrained by various cosmological considerations. 
In case that the Universe is dominated by the saxion, the axion energy density produced by 
the saxion decay should be less than that of one neutrino species to be consistent with the 
Big Bang nucleosynthesis. In addition, LSPs from the saxion decay should not overdose 
the Universe. To satisfy these constraints, one needs to enhance the saxion coupling to the 
SM particles. This is naturally achieved when the jJL term is generated by a superpotential 
interaction between the axion supermultiplet and Higgs doublets. 

This paper is organized as follows. In section 2, we examine how the saxion direction is 
lifted in the presence of mixing between the messengers that transmit SUSY breaking and 
PQ symmetry breaking to the MSSM sector. We then show in section 3 that the model, 
where the PQ scale is radiatively stabilized, can naturally generate the correct mass scale 
not only for ji but also for Bjj, in gauge mediation. Sections 4 and 5 are devoted to the 
discussion of phenomenological and cosmological aspects. We will examine the pattern of 
sparticle masses, the decay of sparticles into axinos or gravitinos, and the cosmological 
constraints on the saxion properties. An important role is played by the saxion/axino 
interactions with the MSSM particles induced by the Higgs coupling to the axion superfield. 
The last section is for the conclusion. 
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2. Axion in gauge mediation 



To invoke the PQ mechanism within the framework of gauge mediation, we introduce heavy 
matter superfields that form vector- hke pairs under the SM gauge groups. These fields are 
classified as 

$ + ^ : SUSY breaking messengers, 
^' + ^ : PQ messengers, 

depending on the way of getting massive. The $ + $ are vector-like also under C/(1)pq 
and directly couple to hidden sector fields that participate in SUSY breaking. They are 
the usual messengers for gauge mediation. On the other hand, the PQ messengers couple 
to the axion superfield S through Yukawa interaction and thus acquire heavy mass after 
PQ symmetry breaking. 

One might think that minimal field content is prepared for the PQ mechanism^ to work 



in gauge mediation. However, previous studies [12, [T^, 14 1 have noticed that there should 
be new interactions transmitting SUSY breaking to the PQ sector in order to fix the PQ 
scale^. This feature is observed under the assumptions that U{1)pq is spontaneously broken 
by a single field, S, and that soft terms only receive gauge-mediated contribution. Hence, 
one may extend the model to include extra SM singlet fields carrying PQ charge, or add an 
additional source of SUSY breaking such as gravity mediation [|l^, Q . Another interesting 
approach we would like to pursue here is to consider the case that some messengers of 
the two sectors, say $ and have the same charge under all the symmetries of the 
theory. Then, there arises mixing between them, which makes S feel SUSY breaking at the 
same loop level as SM-charged scalars do through gauge mediation. Since SUSY breaking 
generates a potential for the saxion, such mixing can play an important role in determining 
the PQ scale. 

In this paper, we consider a simple axion model^ within minimal gauge mediation where 
the messengers <1> + $ and ^' -f- ^ belong to 5 + 5 representation of the SU(5) into which 
the SM gauge groups are embedded. The gauge coupling unification is thus preserved. To 
allow mixing between the 5 messengers, we simply take the PQ charge assignment such 
that ^ carries a charge opposite to that of S while all the other messengers are neutral. 
The model is then described by the superpotential 

W = Wo{X) + y^X^^ + y^S-^)^ + yxX^^ + ysS"^^, (2.1) 



^To solve the strong CP problem, the PQ symmetry should be anomalous under QCD interactions. If 
4' + 4' are charged under QCD, one obtains a KSVZ-type (hadronic) axion model |^|. A DFSZ-type axion 
model ^] is otherwise obtained for the Higgs bilinear HuHd charged under f/(l)pQ. 

^See also for other interesting observations on axions in gauge mediation. 

^Though we are assuming that the model has the PQ symmetry to solve the strong CP problem, it is 
also possible to consider other cases where S corresponds to a flaton field driving thermal inflation. In such 
case, C/(l)pQ needs not be exact, but C/(l)pQ-breaking terms should be small enough so that the potential 
for S can remain approximately flat in the supersymmetric limit. Most of our discussion can apply to these 
flaton models The property of the angular scalar of S would however be quite different because its 

mass is sensitive to the ?7(l)pQ-breaking terms. 
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in the field basis where the Kahler metric is diagonal, ignoring Planck-suppressed operators. 
The above superpotential includes all renormalizable couplings consistent with the SM 
gauge invariance and the PQ symmetry. The effects of hidden sector SUSY breaking are 
parameterized by a background singlet field X, whose Kahler potential should be included 
to correctly compute the anomalous dimension of operators depending on X. 

The PQ symmetry ensures that S corresponds to a flat direction in the supersymmetric 
limit for all the messengers fixed at the origin. Transmitted to the PQ sector by + 
the SUSY breaking effects will lift this flat direction and fix the vacuum expectation value 
(VEV) of the saxion, i.e. the PQ scale. Meanwhile, the axion remains massless until the 
QCD instanton effects are turned on. To examine the saxion potential, it is convenient to 
construct an effective theory with messengers integrated out. Here we assume 
so that the mass of <l>-|-$ is dominated by supersymmetric contribution. Before proceeding 
to the analysis, one should note that the 5 messengers can be redefined further, without loss 
of generality, so that either X^^ or is removed in the superpotential while keeping 

the Kahler metric diagonal. It is thus natural to assume 

'^^'+M = o(l), (2.2) 



in the canonical basis, but and y<j, may be hierarchically different from each other. The 
theory has thresholds at scales A$ = y<i>|X| and Aij, = yiiilSI. For fixed A$, redefining the 
5 messengers appropriately, messenger mixing can be treated perturbatively along the flat 
direction with A,j, far from the scale A$. 

2.1 Saxion potential 

Let us first examine the saxion potential at the region with Aij, = y^\S\ <^ A$. To derive 
the effective action for S, the SUSY breaking messengers can first be integrated out taking 
the field basis where the superpotential is written as 

W = Wo{X)+y^X<^^ + yq,S'^^ + ysS'^^, (2.3) 

for the canonical Kahler potential. This approximately corresponds to the mass basis 
where the term gives small mixing between 5 messengers. The $ + ^, which are 

integrated out below the scale A$, communicate the SUSY breaking to SM-charged fields 
via gauge interactions. In addition, the communication does occur also through the Yukawa 
interaction in two ways: 

1. The anomalous dimensions of S and ^ are discontinuous at A$. As a consequence, 
the saxion acquires soft mass at the two-loop level |l^ 



STT^dyl 



dlnQ 



^ dlnQ 



~ N^N^ [5{{5N^N^ + iV$ + N^)yi + yi) - 2{8gi + 3gi + gi)jyi^, (2.4) 

where we have neglected the splitting between the Yukawa couplings of doublet and 
triplet messengers. The scalar components of ^ -|- ^' also receive additional soft mass 
terms at two loops, besides the ordinary gauge-mediated contributions. 
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2. The effective Kafiler potential for ^ receives a correction 



5K = 



vim 



(2.5) 



from the tree-level exchange of the SUSY breaking messengers. Hence, soft mass 
terms for the scalar components of ^ receive additional contribution 



Mi 



-iV$(167r 



2\2 



A| 



(2.6) 



The above contribution is always tachyonic and becomes sizable or even more domi- 
nant than the gauge- mediated soft mass for a saxion value close to A<j,. 

Here Q denotes the renormalization scale, (iV,j,) is the number of <J? + $ (^ + ^f) 
pairs, and Mq sets the overall scale of soft terms generated by the loops of SUSY breaking 
messengers 



Mn 



1 

167r2 X ' 



(2.7) 



Note that trilinear couplings for the scalar components of ^' + ^' also receive contribution 
mediated through the Yukawa coupling ys- The explicit expressions for soft terms are 
given in the appendix. 

With the knowledge of how SUSY breaking is transmitted to the PQ sector, we further 
integrate out the remaining messengers \1' + ^ under a large background value of S. The 
effective action is then determined by the running wave function of S 



£eff= / d^e ZsiQ = y^\S\)\S\^ 



from which the equation of motion for reads 



F^ 



g -"2(^5 (A*) -7s (A*)) ^ 



(2. 



(2.9) 



neglecting corrections suppressed by A|, / A| . Here are the anomalous dimensions above 
and below A$, respectively. Hence, the scalar potential for the saxion is generated as 



y = Fo + m|(Q = y*|5|)|5p 



(2.10) 



where m| is the running soft mass of S in the theory between Avp and A$, and a constant 
Vo has been added to cancel the cosmological constant. It should be noted that rn^^Q) 
depends on |5| itself because the Kahler correction ( p.5| ) generates soft terms for the scalar 
components of ^' + ^' that affect the running of ?7i|. Thus, in the region with y^ilS*! <^ A$, 
the saxion potential has a slope approximately given by 



dV 



2Mi\S\ d\S\ 



m^(A^) ^ 5N^N^ 



Mi 



87r2 



2y| 



levr- 



\ y^\s\ 

Aff> 



In 



m>\s[ 

A* 



(2.11) 
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where the Yukawa coupUngs are evaluated at A$, and C^r depends on y| g and gauge cou- 
plings. The logarithmic dependence originates from running between A<j, and Aq, through 
the Yukawa interaction with ^' + Due to the radiative effects from the color-charged 
scalars of ^' + Cq, has a positive value of order unity. 

The slope of potential ( p. 11 ) shows that the saxion can be stabilized by the balance 



between two effects, i.e. the SUSY breaking mediated by $ + ^ at A$, and the renormal- 
ization effect through the Yukawa interaction with \I' + The second part of the slope 
monotonically increases as a function of |5| for y*|5| < O(0.1A$), and crosses zero at 
yitlS'l = 0(A$/87r2) because the two contributions in the bracket have the opposite sign. 
Hence, depending on the value of m| at A$, the potential develops a minimum along the 
saxion direction as follows. If m|.(Aj) has a positive value of 0{Mq) or less, the saxion is 
stabilized at y^\S\ < ©(A^/Svr^). On the other hand, for negative m|(Aj), a minimum 
appears at a scale rather close to A<j, where the Kahler correction ( |2.5D becomes important. 
In fixing the VEV of the saxion, the crucial role is played by the messenger mixing as can 
be seen from that the potential only has a negative slope in the limit — )■ 0. 

Let's move on to the opposite region with A^ ^ A$ along the saxion direction. In 
this region, the correct procedure for constructing the effective theory is to first integrate 
out the PQ messengers at the scale Aij,. For this, we take the field basis such that 

W = Wo{X) + y^X^^ + y^S"^^ + yxX^^, (2.12) 

for the canonically normalized fields. Integrating out \I' + ^', one obtains the effective action 
for X determined by its running wave function. In the effective theory, S does not have 
any renormalizable interactions, and the equation of motion for gives 

pS 2 |j^|2 

— ~ -57V$7Vv,^^^Mo, (2.13) 

where corrections suppressed by A|,/A|, have been neglected. Hence, the leading contribu- 
tion to the saxion potential comes from the dependence on S of the effective wave function 
of X. Because the anomalous dimension of X is discontinuous at the scale A^, the slope 
of the potential is derived as 



dV 5N^N^ 



2 _C<i> 2 



yx - T^Zr^y^ 



F^|2, (2.14) 



where a positive constant C<j> = 0{g'^) parameterizes the contribution induced at the three- 
loop level. The potential thus increases at y^\S\ ^ A$ as a function of l^l, unless yx is 
smaller than 0{y^/8Tr'^). This property is cosmologically favorable because the saxion may 
be displaced far from the minimum at the end of inflation. If this happens, the positive 
slope will make the saxion roll down toward the true minimum. 

The relation ( ^.14 ) also gives information about the potential at saxion values close 



to A$, for which messenger mixing can no longer be treated as a perturbation. Instead of 
constructing an effective theory, we use the property that the slope at A^ » A$ is positive 
for yx = 0{y^). This implies that there must exist a minimum below or near A$ since 
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m| is driven negative at low scales by the Yukawa coupling with ^' + For the saxion 
stabilized far below A$, the vacuum structure can easily be examined treating the mixing 
perturbatively. Note that the potential at y^\S\ <C A$ is essentially determined by yx^/^s at 
A$ and insensitive to the details of the potential at large Another possibility is that 
a minimum lies close to A$, which generically requires a rather small y|. The existence of 
minimum is ensured by the positive slope at A<j, » A$. 

It is worth discussing the situation that there is no mixing between the messengers, 
as usually assumed in gauge mediation. This corresponds to the limit that yx,s vanish. 
Prom the relations ( |2.1l| ) and ( |2.14| ), one then finds that the potential runs off to infinity 
along the saxion direction. Hence, additional SUSY breaking effects are needed to stabilize 
the saxion. A natural candidate for this is gravity mediation since the saxion potential is 
generated at three-loop level. Indeed, a higher dimensional operator oc |X|^|5p/Mp; in 
the Kahler potential gives the gravity-mediated contribution 



6V 



(2.15) 



which can compete with the gauge-mediated one to stabilize the saxion for a positive k of 
order unity jl3| ]. In the presence of messenger mixing, however, the saxion potential has 
a slope as ( p. 11 ) at \S\ <C A$ as long as X is smaller than O{10~^Mpi), and the above 
contribution becomes important only at |5| > 0{yxMpi/V8iT^k). 

We complete this subsection by summarizing the role of mixing between SUSY breaking 
and PQ messengers. Such mixing indicates that there exist some SM-charged heavy fields 
that directly couple both to the SUSY breaking fields and to the PQ breaking field. It is 
the SUSY breaking effects transmitted by these fields that radiatively generate a potential 
for the saxion and fix the PQ scale. Moreover, the mixing prevents a runaway behavior of 
potential at large saxion values in gauge mediation. This would be cosmologically relevant 
for the saxion to settle down to the true vacuum. 



2.2 Vacuum structure 

Messenger mixing can be treated as a perturbation at y^\S\ <C A$ to construct the effective 
theory for S. Taking into account that ^ receives a Kahler correction ( |2.5D that contributes 
to its effective wave function, we examine the vacuum structure focusing on the case that 
the saxion is stabilized at yiiilS"] < ©(A^/VSvr^). From the effective action for S, the saxion 
a and the axino a are found to acquire SUSY breaking mass as 



Mo 



47r2 



C^yl 



8^ 



42/1 



IGvr- 



.y^SQ 
Aff, 



In 



A* 



+ 



IGvr- 



./y^So 
A$ 



, / y^So 



(2.16) 



for the axion superfield expanded around its VEV, 5 = {So+a/\/2)e''''/^'^'' + V29a+e'^F^ . 
Here the couplings are evaluated at Q = yqiSp, and Aq^i = O {10 N^N^ygMp) are the 
trilinear couplings associated with the Yukawa couplings of the PQ triplet and doublet 
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Log,„|S| (GeV) v,,,(Aj) 

Figure 1: Saxion stabilization in the model with N^, = 1, 7V$ = 2, A$ = 10^'^GeV and ?/$(A$) = 
0.3. The left panel shows the scalar potential along the saxion direction: the solid line is for the 
case with (y*,?/s) — (0.3,0.5) at A$, while the dashed one is for (y*,?/s) — (0.5,0.4). In the 
right panel, one can see the dependence of the saxion VEV on the Yukawa couplings y^^g at A$: 
= 109-i°GeV (red), = 10i''-"GeV (green). So = lO^-i^GeV (blue). The marker + means 
m|(Aj) > 0, while — for the case with m|(A^) < 0. 



messengers, respectively. For 2/|, 5 = 0(0.1), which are the plausible values, the saxion 
mass lies in the range 

O <m.<0 (^,/M8^Mo) , (2.17) 

where the upper bound is obtained when i/^Sq = ©(A^/VStt^). Though radiatively stabi- 
lized, a can be as heavy as the color-charged MSSM sparticles. This is because messenger 
mixing induces a correction to the Kahler potential for ^ whose loops contribute to the 
saxion potential. On the other hand, the axino mass is rather insensitive to the Kahler 
correction ( |2.5| ) , and has the value 

m, = (^) , (2.18) 

for y"^ g = 0(0.1), thereby lighter than the MSSM sparticles. 

To illustrate how the PQ scale is fixed after SUSY breaking, we provide some examples. 
Fig m shows the scalar potential along the saxion direction in the model with Nq, = 1, 

= 2, A$ = lO^^GeV and y$(A$) = 0.3. Depending on the Yukawa couplings yq,^s at 
the scale A$, the saxion is stabilized in the following ways: 

= (0.3,0.5) : So ~ 1.3 x lO^OCeV, Vq ~ {OAMofSl ~ O.SMq, 
{y^,ys) = (0.5,0.4) : So ~ 4.2 x lOi^GeV, Vo ~ {0.9MofSl m„ ~ I.TMq, (2.19) 

where m|(Aj) is positive for {ym,ys) = (0.3,0.5), while negative for the other case. For 
m|(Aj) < 0, the saxion potential develops a minimum at a scale near A$ while providing 
a rather large mass to the saxion. As discussed already, the Kahler correction ( p. 5]) to the 
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PQ messenger generates a potential for the saxion that becomes important at saxion values 
close to A$. In Fig |l[ one can also see how the saxion VEV is fixed depending on y^/^s at 
A$. For ys < 0.42 where m|.(Aj) is negative, the contribution from the Kahler correction 
(12. 51) stabilizes the saxion at |5| > ©(lO-^A^). It is also possible to obtain Sq = ^^(A^), 
for instance, by taking y^, = 0.2 and ys smaller than 0.3. In this case, the eff'ective theory 
for ^' + ^ constructed by integrating out <I> + ^ is not reliable. Nonetheless, the relation 



(2.14) tells that the potential has a minimum as long as yx = 0{y^). 



3. Higgs /i and Bfx terms 

Although it is an attractive mechanism for generating flavor and CP conserving soft terms, 
gauge mediation requires some additional structure to account for the origin of the fi and 
Bfj, terms in the MSSM. In particular, it is quite non-trivial to obtain an acceptable value of 
B in theories with gauge mediation. If one introduces a direct coupling of the Higgs bilinear 
HuH^ to the SUSY breaking field X in the superpotential, one obtains B = 0{8tt'^ AIq) and 
thus needs an unnatural fine-tuning to achieve the electroweak symmetry breaking. One 
may instead consider an effective Higgs coupling in the Kahler potential 



I d^ef{X,X*)H^Hd, (3.1) 



which relates ^ to the SUSY breaking parameters^. However, this operator generically 
gives B of ©(Stt^Mq) again. To avoid large B, f should have a particular dependence on 
X such that generates fi but not Bfi [19|, unless there are other SUSY breaking fields. For 



example, one can use / = X* /A with some mass scale A |2C|. The dynamics that connects 
the Higgs sector to the SUSY breaking sector in such a particular way would generally 
affect other MSSM soft terms generated by gauge mediation. 

Here, we take an alternative approach to solving both the ^ and Bfj, problems, which is 
provided by the PQ mechanism incorporated into the gauge mediation. In fact, a natural 
solution is to consider a coupling between HuHd and the axion superfield S. The /i term 
is then induced only after the PQ symmetry is broken. Furthermore, because the radiative 
stabilization of the saxion leads to 

pS 

— ~ -5iV$iV^y|Mo, (3.2) 

Do 

a generic coupling of HuHd with S is naturally expected to give B of the order of MSSM 
sparticle masses for y| = 0(0. 1)^. The above relation is a consequence of mixing between 
SUSY breaking and PQ messengers. In the absence of such mixing, though the saxion can 
still be stabilized by adding additional SUSY breaking effects such as gravity mediation, 
/So would have a value much smaller than Mq since S couples to X at more than 
two-loop level. The relation (3.2) implies that there are simple mechanisms operative to 



generate fi and Bfi terms required for proper electroweak symmetry breaking: 



*This is a generalization of the Giudice-Masiero mechanism jl^] that generates fj, by SUSY breaking 
eflect. 

similar idea to suppress B was considered by [pi| in a different model. 
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Kim-Nilles (KN) mechanism^ : After the PQ symmetry is broken, the non-renormahzable 
term in the superpotential 



02 



/:kn = J d'dXj^HuHd + h.c. (3.3) 
generates ^ and terms with 

For the saxion stabihzed at a scale around lO^^GeV, which is well within the invisible 
axion window consistent with astrophysical and cosmological bounds, ^ is generated 
at the soft mass scale with A = 0(0.1). In addition, we obtain B of the correct order 
of magnitude. 

• Giudice-Masiero (GM) mechanism : For the Higgs fields that couple to S through 
the effective Kahler potential term 

Cgm = J d^On^HuHd + h.c, (3.5) 
both /i and B are induced by SUSY breaking effects 

^ = K—, ^ = ^- (3-6) 

The /i and Bfj, terms are thus of the desired size. The above coupling in the Kahler 
potential can arise, for instance, by integrating out heavy fields in the model with 
the superpotential terms TiHuHd + TiSS' for S' having a wave-function mixing with 
S in the Kahler potential. 

It is important to note that, since the phase of /Sq is aligned with Mq, the Bfi term 
does not introduce new source of CP violation in either mechanism. Notice also that the 
Higgs coupling to S fixes the PQ charges of the MSSM matter fields. In order for a to 
play the role of the QCD axion, C/(1)pq should be anomalous under QCD. This requires 
non-zero value for the QCD anomaly coefficient, N = N^i, it 6, + for the GM while — for 
the KN mechanism"^. 



®A non-renormalizable superpotential coupling of H^Hd to the PQ breaking field has been considered to 
explain the /i term in various SUSY breaking schemes . The size of B however depends on the mechanism 
stabilizing the PQ scale. 

^Models that incorporate the KN mechanism with A'^ — ±i are free from the domain wall problem. For 
example, in the model with Ni^j = 5 and Nif, — 1, one obtains TV = — 1, and the MSSM gauge couplings 
remain perturbative up to the unification scale as long as fO^GeV < A* < 0.fA$. Meanwhile, for the 
case that /i is generated by the KN or GM mechanism with A'^ 7^ ±i, one can consider the situation that 
the saxion is displaced far from the origin during the inflation. Then, the reheating would not restore 
U{l)pQ after inflation. In addition, provided that the fluctuation around the initial displacement due to 
the quantum fluctuation during the de-Sitter expansion is small enough, the saxion will settle down to one 
of the I A'' I degenerate vacua. This will provide a solution to the domain wall problem (see jiij, for a similar 
consideration) . 
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Another important consequence of the Higgs couphng to S is that mixing between 
the saxion (axino) and neutral Higgs (Higgsino) fields is induced after electroweak symme- 
try breaking. Through the mixing, the saxion and axino interact also with other MSSM 
particles. The saxion mixing term with neutral Higgs bosons is obtained from 



by making the replacement 



1^1 



5*0 V2 



with Ca 



din |/i| 



dln\S\ 



(3.7) 



(3.8) 



S=So 



where we have used that / S does not depend on S, as can be seen from (2.9). Therefore, 
the saxion slightly mixes with the neutral CP even Higgs bosons through the interaction 
suppressed hy v/Sq with v'^ = {\H^\)'^ + {\H^\)'^. Similarly, the axino has tiny mixing with 
the neutral Higgsinos determined by 



(i7°iJ,° + i?Ma + h.c 



uO\~ 



Sn 



with 



C-a 



(3.9) 



(3.10) 



S=So 



dlnS 

Note that the coefficient Co- crucially depends on the mechanism for generating the n term: 



KN 



c, 



a\GM 



Q- 



dQ 



Q=S/*S'o 



O 



V87r2 



whereas the size of Ca is insensitive to the form of the Higgs coupling to S: 

CalxN = 2, CqIgM = —1- 



(3.11) 



(3.12) 



Let us examine the saxion/axino couplings to the MSSM sector. First, there are the 
interactions, (T|i/^| |-ff^ ^| + aH^H^ + ^H^Ji, that are induced directly from the Higgs 
coupling to S. The couplings for these interactions are non-vanishing even in the limit 
u — )■ 0. To derive other interactions between the saxion/axino and the MSSM particles, the 



small mixing terms from (3/7) and ( p.9[) should be removed by performing an appropriate 
field redefinition. Indeed, the saxion and axino couplings can be read off from the MSSM 
Lagrangian by the substitutions 



CaV 1^1 



d,u 



So ml 

CgV 

So 



mi: 



V2 ' 



(3.13) 
(3.14) 



where mh is the mass of the lightest CP even neutral Higgs boson h, and A'^^ are non- 
vanishing because ^ mix with the bino and neutral wino. The mixing parameters are 
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presented in the appendix. After the replacement, one must diagonahze the mass matrix 
for {B,W^ , H^, H^). In addition to those proportional to Ca^a, interactions between the 
saxion/axino and MSSM particles are also generated by the loops involving PQ messengers 
that acquire mass from the VEV of S. These couplings can be derived from the dependence 
on S of the MSSM gauge couplings 

-5- = ^ln(S'*5) + (5-independent part), (3.15) 

in the effective theory with ^ + 'I' integrated out. The saxion decay to MSSM particles 
and the decay of heavy sparticles into axino will be discussed later, after examining MSSM 
sparticle masses. 

We stress that the axion superfield can play an important role not only in solving the 
strong CP problem but also in explaining the presence of the // and Bfj, terms within theories 
with gauge mediation. This nice feature stems from the mixing between SUSY breaking 
and PQ messengers, through which the SUSY breaking is communicated to the PQ sector 
and radiatively stabilizes the PQ scale. In fact, an appropriate PQ charge assignment is 
the only thing that was needed to allow such mixing between the messengers. 



4. Phenomenological implications 

In this section, we discuss the phenomenological aspects of the model. 
4.1 Sparticle masses 

To derive soft terms for the MSSM fields, we begin by summarizing the possible range of 
threshold scales in the theory. The PQ scale, which is radiatively stabilized in the presence 
of mixing between messengers, is constrained by various astrophysical and cosmological 
observations ||5|. On the other hand, an upper bound is put on A$ = ?/<i>|A| to suppress 
gravity mediation that in general generates flavor- violating soft terms with size F'^ /Mpi. 
These constraints lead to 

lO^GeV <So< lO^^GeV, lO^o < \X\ < lO^^GeV, (4.1) 

where the lower bound on X has been put to concentrate on the case that messenger mixing 
can be treated perturbatively. Since the theory contains heavy messengers, the MSSM soft 
terms at TeV scale are determined by the parameters {Mq, A$, iV$, Aij^, A^^}, while n and 
B are generated by the KN or GM mechanism. 

At the higher threshold A$, the SUSY breaking is transmitted to the MSSM sector 
by <1> + ^> through ordinary gauge mediation Q. The threshold effects induced by these 
messengers generate gaugino and scalar soft masses as 

4^ = 2A$Cf5.'(A^), (4.2) 
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for being the quadratic Casimir of the corresponding field. Soft trihnear terms arise 
at two- loop order, and thus are negligible at A<j,. Since the radiative corrections due to 
the SM gauge interaction play the dominant role for the mediation, soft terms preserve 
flavor and CP. The gauge-mediated soft parameters are subsequently RG evolved down to 
= y^So in the presence of SM-charged fields ^' + ^. 

Low energy soft terms below Aij, can be computed integrating out ^ + ^. Since the 
saxion IS stabilized with F^/Sq = 0{Mq), the PQ messengers induce negligible threshold 
effects for gaugino masses and trilinear couplings. However, depending on the saxion VEV, 
soft mass terms for the MSSM scalar fields can receive sizable threshold corrections. This 



is because the PQ scalars acquire soft mass also from the Kahler correction (2^). The 
threshold effect provides flavor-universal soft scalar masses 

^.4A.,^.[C».^-r..|.H]4|lln(^). ,4.3, 

where the gauge couplings are evaluated at Q = A^, , and corrections suppressed by the two- 
loop factor have been neglected. The above contribution includes the hypercharge trace 
term, which is the part proportional to the hypercharge YJ. This term is non- vanishing due 
to the Yukawa splitting of doublet and triplet components in ^ + 

y|(A<,) A$ ^^-^^ 

with ?/| = + y|^)/2, and Mgut being the unification scale. The other term in the 
bracket arises from the loops of ^ + ^ since the supertrace of their mass matrix is non- 



vanishing |2^. For example, depending on the values of y-^^s at Mqut, the model with 
N^ = N^, = 1 and A$ = lO^^GeV leads to 

{y^,ys) = (0.34,0.68) : i? ~ 14.5, Sq ~ 3.5 x lO^GeV, 

{y^,,ys) = (0.24,0.48) : R ~ 14.3, - 8.7 x 10^°GeV, (4.5) 

where we have assumed ^ 1 to evaluate the running of y^^^s from Mqut to A$. In the 
case with {yqi,ys) = (0.24, 0.48) at MquTj the PQ scale is close to A$, and thus the MSSM 
scalars receive non-negligible threshold correction Am? ~ — 0.131^ Mq at A^ ~ 2xlO^^GeV. 
It is worth noting that one can simply change the sign of R by considering mixing between 
5 messengers instead of between 5 ones. Hence, R can be of either sign. 

The threshold effect at for soft scalar masses becomes important when the saxion 
VEV is close to A$ , as can be deduced from its origin. For instance, ify^So = 0{A^/V Svr^), 
the hypercharge trace term makes Am? comparable to the gauge-mediated soft mass for 
scalars that are charged only under U{1)y- This indicates that the lightest ordinary spar- 
ticle (LOSP) can be provided by the stau even for small N^, if R is positive and the saxion 
is stabilized at a scale near A$. However, for the saxion VEV much lower than A$, the PQ 
messengers only give negligible threshold to soft terms. Meanwhile, the gauge-mediated 
sfermion masses satisfy two sum-rules 



^ Ym^i =Y.{Bi- Li)mj = 0, (4.6) 
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where the sum is over one generation of sfermions, and B (L) is the baryon (lepton) number. 
Affected by the PQ threshold Am?, these sum rules can be used to extract information 
about the PQ sector. 

Finally, we examine which particle is the LSP. The theory contains two light fermionic 
sparticles, the axino and the gravitino. To determine the gravitino mass, we assume that X 
is the goldstino superfield whose F-term cancels the cosmological constant. The gravitino 
then absorbs the fermionic component of X, X, to become massive with 

Mo ^ Mpi' ^ ■ ^ 

which lies in the range 10"'^ Mo < m3/2 < 10~^Mo for the theory with (j^. Using that 
the axino acquires mass of ©(Mo/Svr^) for y|, 5 = 0(0.1), one also arrives at the relation 

"^3/2 X 

m-a ^ IQi^GeV ^ ^ 

This shows that the axino becomes the LSP if X > O(10^^)GeV. The LSP would otherwise 
be given by the gravitino. Though we will not consider it here, there is a possibility to have 
an axino LSP even for X < lO^^GeV. One way is to consider y'^ ^ <^ 0.1. For instance, 
in models with y| = 0(10"^)^ obtain m^ = ©(IQ-^Mq) and m^ > O(O.lMo). For 
y| ^ 0.1, a Higgs coupling to S will give \B\ <C 0{Mo), with which it is still possible to 
achieve the correct electroweak symmetry breaking. One can also consider other models 
where X is not the main component of the goldstino. Then, F-^ will have a VEV less than 
0{m^l2Mpi), and thus the gravitino can be heavier than the axino for X < lO^^GeV. In 
this case, one would need 171^/2 ^ "^o- iii order not to destabilize the PQ scale, because 
gravity mediation provides soft mass typically of O {171^/2) to the saxion as well as to the 
MSSM scalars. 

4.2 Decay of sparticles 

The heavy MSSM sparticles rapidly decay into the LOSP, denoted by x, which subsequently 
decays into lighter sparticles, i.e. into the axino or gravitino. The decay of x occurs more 
slowly because the axino and gravitino are very weakly coupled to other particles. Here we 
are assuming R-parity conservation. Measurement of the decay length of x will give direct 
information either about the SUSY breaking scale or about the PQ scale, depending on 
which of axino and gravitino is the main decay product^. In fact, the axino and gravitino 
have similar type of interactions that mediate the decay of heavy sparticles. 

At energy scales much higher than m3/2, the gravitino G effectively behaves as a 
goldstino. To study the decay of heavy sparticles into gravitino, the effective interaction 



Lagrangian for the goldstino component can be written in non-derivative form |25|. Using 
the relation F^ = —IQir"^ MqX ~ y/j>m^j2Mpi, one obtains 

- ' ( '^^^~'^K *ilX + -^^Xa'^-XF ]+hc (4 9) 

16^2^ I Mo '^^^ + 4^/2Mo '^'^i^ ' ^ ^ ^ 



^However, the decay length alone would not allow us to distinguish between SUSY breaking scenarios 
where the LSP is given either by the axino or by the gravitino |24[. 



- 14 - 



independently of how the SUSY breaking is mediated. Here A stands for the gaugino, and 
F^i, is the corresponding field strength, while <j) is the scalar, and is its fermionic partner. 
The gravitino interactions are proportional to the mass splitting in the supermultiplet, and 
inversely proportional to . 

Since the // term is generated after the PQ symmetry breaking, Ca has a value of order 
unity, and the axino interacts with MSSM particles through 

KMu-a + {KdBHld + ~rfKd + ~rBf + rw'f + H%a^Hl,Z,) , (4.10) 

where we have omitted the couplings. For those in parenthesis, the axino couplings are 
obtained after the replacement ( |3.14 ), and thus would vanish in the limit u — )• 0. Here 
is the Z-boson with mass Mz, while / denotes the SM fermion, and / is its scalar partner. 
In addition, there are axino interactions induced at the loop level 



(4.11) 



whose couplings are determined by the dependence on S of the gauge couplings after 
integrating out ^' + 'I'. Because of the axino mixing with ^, the above couplings also 
receive contribution from the loops involving MSSM particles that become massive after 
the electroweak symmetry breaking. 

Using the Lagrangian for axino/gravitino interactions, one can estimate the lifetime 
of X) which is also subject to cosmological constraints. In the model, the LOSP can be 
provided by the stau or bino. Let us first examine the case that the bino is the LOSP. The 
coupling for the interaction hBa receives contribution from that of H^*^BH^ ^ due to the 
axino component of neutral Higgsinos, and from d^dvP' through the mixing between B 
and H^^. This coupling is of 0(Mz / Sq), and thus the bino mainly decays into h and a 
with decay width 



1 f \ / IQiQGeV y 

10-8sec V200GeVy \ Sq ) ' 



(4.12) 



for Mj^ > nifi + rria = rrih + ©(Mq/Stt^), assuming that other neutral Higgs bosons are 
very heavy. It is easy to see that the bino decay into gravitino is highly suppressed. For a 
bino with < nih + rrid, the dominant decay channel is i? — >• Zd. The coupling for this 
decay process is additionally suppressed by 0{Mz/n) compared to that of hBd, because it 
requires bino-Higgsino mixing as well as the axino-Higgsino mixing. As it is mediated by 
the interaction induced at the loop level, the decay S — t- 07 has a small branching ratio. 
On the other hand, in the case of a stau LOSP, the decay takes place via the interaction 
f*ra with coupling of 0{{MzN'l + rrifN'-^Q/ cos/3)/5o). The decay rate is estimated as 



10- 



' sec 



Mz cos 2/3 
0.3M^ 



(l-natan/3)- 

102 



/ ruf \ /lOi^GeVV -,o^ 
(200Giv)(^^J ' (^-^^^ 



with n-a = 0{Ml/ixM^), and tan/3 = Therefore, for a bino or stau LOSP, 

the LOSP will decay mainly into axinos, rather than into gravitinos. Measuring its decay 
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length, one can thus extract information about the PQ scale. It is interesting to note 
that, depending on the PQ scale, a bino LOSP can decay inside the detector while leaving 
displaced vertices. For a stau LOSP, there is a possibility that decaying staus can appear 
in the detector often enough to measure their charged tracks. Since x ii^s a lifetime much 
shorter than a second, the nucleosynthesis does not place any significant bound. 

The LSP, into which all the sparticles eventually decay, can be either the axino or the 
gravitino depending on the scale of gauge mediation, and constitutes the dark matter of 
the Universe. If it is lighter than the gravitino, the axino becomes a good candidate for 



the cold dark matter [26, 27 1. The axino can become the LSP when X > O(10^^)GeV for 
5 ~ 0(0.1). In this case, the gravitino decay occurs through the interaction 

Ant = ^^^^^^^'^^^d.a + h.c, (4.14) 

which leads to 

r 1 ( Sn'^myX f Mq y 

G^aa ~ 5 ^ ^QiSggg ^ Mo J VsOOGeV/ ' ^' ' 

The gravitino decay will thus produce an axino and axion. On the other hand, if heavier 
than the gravitino, axinos produced by the decay of x will decay into gravitino. Since it is 
mediated by the effective interaction 



^nt = ^^^i^^^^'^^M" + h-c, (4.16) 
for rria > rn^i2-, the decay a — )• Ga occurs with 



r 1 / ma y/ SA^ y/ Mo ^ 

^^^•^ 0.5 X 10i3sec Vm3/2y \ ) VsOOGeVy ' ^' ' 



where we have used the relation ( |4.7| ). The axino will decay with a long lifetime, producing 
gravitinos together with axions. Note that late decay of axino/gravitino produces LSPs and 
axions, which may be warm or even hot at present unless a and G are highly degenerate in 
mass. In fact, having a free-streaming length much larger than O(10)Mpc, LSPs produced 
by such late decays will behave like a hot dark matter. The energy density of hot dark 



matter is severely constrained by the CMBR and structure formation [28, |29[]. We will 
return to this issue in the next section. 



5. Cosmological aspects 

The theory contains the saxion that has a rather flat potential generated after SUSY 
breaking and interacts with other particles with coupling suppressed by the PQ scale. This 
scalar may play some non-trivial role in cosmology as its potential can receive additional 
sizable contribution at early Universe. The relic abundance of dark matter depends on the 
cosmological evolution of the saxion. It is thus of importance to understand the saxion 
properties. 
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5.1 Saxion decay 

Because it acquires mass as 0(0. IMq) < rria < 0{Mq) depending on S'o/A$, the saxion 
will decay into axino, gravitino, and light MSSM particles. The interaction relevant for its 
decay can be derived from the effective action (|2.8| ) 



"int 



(5.1) 



where we have included the effective interaction with the goldstino. The coefficient — 
{dlnma/dln\S\)\s=So has a value of 0(0.1) or less for i/'^Sq < O(10~^A$). If the saxion 
is stabilized at y\ii\S\ = O(0.lA$), the Kahler correction ( |2.5| ) becomes important and 
leads to Ca = 0(1) and mo- = 0(Mo). Besides the above interactions, the saxion also has 
couplings to MSSM particles when the fi term arises after PQ symmetry breaking: 

a^o^o + a\H'jHl,\ + [hIJ f + hI*,BH% + H^Jf + H^Z^^Z^ + •••), (5.2) 

where the first two terms come from the Higgs coupling to S, while the other interactions 
require the saxion component of ^. The ellipsis contains the interactions with Wj^ and 
W^'^. There are also saxion interactions induced at the loop level, which include 

aJr+aFf^'^F^,, (5.3) 

where the couplings can be derived using the relation (|3.15|), and making the replacement 



(3.13). Though it is suppressed by the loop factor, the interaction aff* can be important 
because the saxion coupling from is suppressed by CaAf /Mq = 0{C(j/S'k'^) in gauge 

mediation. Here ^/ is the soft trilinear coupling. 

Before examining the partial decay widths of the saxion, we note that the color-charged 
sparticles are much heavier than uncolored sparticles in minimal gauge mediation. This 
results in that the electroweak symmetry breaking is achieved under a fine tuning of a few 
percent between the Higgs mass parameters. So we shall consider the term with^ 

M2 ' M2 ' M2 ^ M2 - ^^-^^ 
For later discussion, it is convenient to define the total decay width of the saxion as 

1 

with Ba being the branching ratio into axions. Let us now examine the saxion decay. From 



(5.1), the decay widths for a — )■ aa and a — )• aG are determined by 



^Using that the combination Ma/gi is invariant under the RG evolution at the one-loop level, one obtains 
Mg ~ Q.22N^Mo at Q = ITeV. The bino mass is thus about lOOGeV for iVj, = 1 and A'/o = 500GeV. In 



the subsequent discussion, we will treat Mg as a free parameter satisfying M? ^ Mq 
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The saxion also decays into MSSM particles via the interactions (5^) and (|5.3| ), and thus 
crucially depends on the form of the Higgs coupling to S. In particular, for = the 
saxion decay occurs with 

2 1^1"^ ( 2 Nf^rn^ml k^M^ \ 

{T^^hh,'^c7^ff,'^a^zz) ^ 4:C^BaT^— k^, — ^ —, — ^ — , (5.7) 

\ (m^ — m^)^ [mj^ — m^yj 

if the processes are kinematically allowed, where kh = 1 + \B\ sin2a/2|//|, kz = sin /3 + 
cos /3, and we have neglected the mass of the decay products. In the decay width into 
the SM fermions, A''^ = 3(1) for quarks (leptons), and kf = N^/ sin/3 for up-type quarks 
while kf = / cos (3 for down-type quarks and leptons. Here k^ j z have a value of order 
unity. Furthermore, if it is heavy enough, the saxion will decay into MSSM sparticles with 
the decay widths roughly estimated as 

where i? is a bino-dominant neutralino with small Higgsino component, and the stau has 
kf = 4iV|C?(5^(Avi/) - gl{MQ))/ha. To estimate the decay width into MSSM sparticles, 
we have used that the interaction aBB arises from aH^H^ and H^*^BH^ ^ due to the 
bino-Higgsino mixing, and that gauge mediation gives A-r = ©(Mq/Stt^). 

5.2 Cosmological constraints 

During the inflationary epoch, the saxion is displaced from the true vacuum because it 
obtains a Hubble-induced mass term. After inflation, as the Universe is reheated, the PQ 
messengers generate a thermal potential for the saxion, SV ~ y|,r^|5p at IS*! <C T. For 
IS*! ^> T, the messengers become massive and only give small thermal effect. Hence, if it 
sits around the origin just after inflation, the saxion is thermally trapped at the origin until 
the temperature drops down to T = 0{Mq). This implies that thermal inflation [|l6| occurs 
when the potential energy Vq ~ m-a^o dominates the Universe. After thermal inflation, 
the saxion begins to oscillate about the true minimum with an amplitude of 0(5*0). It is 
also possible that the saxion is shifted far from the origin during primordial inflation. In 
this case, the coherent oscillation of the saxion starts with an amplitude less than Mpi 
when the Hubble parameter becomes comparable to Since it behaves like matter, the 
oscillation energy would dominate the energy of the Universe if the initial amplitude is 
large enough. 

Let us examine cosmological constraints of the saxion properties in the case that the 
saxion dominates the Universe at an early time, which is a plausible possibility as discussed 
above. Since a late-time entropy production would alleviate the constraints placed on the 
saxion, we further assume that there is no additional entropy generation after the saxion 
decay. An important constraint then comes from the axions produced by the saxion decay 



as they behave like neutrinos [30|. If there is no coupling between HuHd and S, the saxion 
mainly decays into two axions with Ba — 1- Obviously, the produced axions would spoil 
the Big Bang nucleosynthesis. However, the situation can be different if the term is 
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generated after C/(1)pq breaking. From (^), one sees that the saxion couplings to the 
MSSM particles can easily suppress Ba below 0.1 for 

ml < (5.9) 

when the /i term is generated by the KN mechanism. We assume that this is the case. Using 
(|5^), one then finds that Ba has a value between 0(10""^) and 0(10"^) for > 2Mz- 
Depending on rria, the decay is dominated by a — )• hh, a — > ZZ, or cr — )• it with t being 
the top quark. On the other hand, if ma < 2Mz, the saxion mainly decays into bottom 
quarks, while giving 0(10^'^) < Ba < 0.1 for < 0.1|/xp. Therefore, one can naturally 
avoid production of too many axions from the saxion decay^*^ . 

There is also a constraint from the dark matter abundance. First of all, the total 
abundance of LSPs should not exceed the measured abundance of the dark matter^^. In 
addition, an attention should be paid to the LSPs produced from late decays of the ax- 
ino/gravitino because they will contribute to the energy density of hot dark matter. The 
current energy density of cold dark matter is r^cDM — 0.2, while the density of hot dark 
matter is bounded from above as r^HDM ^ 10^'^ |3^, In the situation under consid- 
eration, the dark matter abundance essentially depends on the decay temperature of the 
saxion 

^ / 90 , ^../0.01\^/V \3/2 /10iiGeV\ , 

for g^{Ta) = O(102), with being the effective degrees of freedom of the radiation. Note 
that there are two main processes for the production of axinos and gravitinos, {i) the decay 
of the saxion cr, and {ii) the decay of the LOSP x- If > 2m^, LOSPs will be produced 
directly from the saxion decay. In addition, there are thermally generated LOSPs when 
is higher than Tf. Here Tf ~ m^/20 denotes the freeze-out temperature of below which 
X decouples from the thermal bath. 

We first concentrate on the direct production of axinos and gravitinos from the saxion 
decay. The yield of axinos produced from the saxion decay is the sum of the axino yield 
from a — >■ aa, Yq and that from a — )• aG, y^i2- 

YI = YS^Y^I2- (5-11) 
Using the relation (5.6) and m-a = ©(Mq/Stt^), each term is evaluated as 



-^(7 3 TJj F^— >.aa 

° ~ 2^ F^ 



A ( BaY'^ ( m~a ^ 2 ^ io2GeV ^ lO^^ GcV ^ .^^ 



^"Models in |3^ have used this property to suppress the axion production from the saxion decay. In the 
model of where the yi term is generated by the GM mechanism, the authors introduced additional SM 
singlet to achieve Ca = 0(1). On the other hand, in axionic mirage mediation which also incorporates 
the GM mechanism, the large entropy released by the modulus decay dilutes the axions produced by the 
saxion decay with Ba — 1. 

The axion also contributes to the energy density of dark matter. Via ~ 0.4e'^(So/10^^GeV)^-^** with Ba 
being the initial misalignment angle. 
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and 

the latter of which is the same as the gravitino yield from a — )• dG. Note that, when 
the saxion is stabilized with y'^S^ < O(10~^A|), one generically obtains Ca = 0(0.1) and 
rrifj = 0(0. IMo). For yij-So = 0(0. 1A$), Ca is of order unity, but the saxion acquires mass 
of 0(Mq). As discussed already, Ba < 0.1 can be achieved for <^ To be consistent 
with the cosmological observations, y^^^2 should satisfy the constraints 

1- < 3.6 X 10-1° (i^) and < 1.8 x 10-- (i^) , (5.14) 

for the axino LSP, where the latter one comes from the constraint on the hot component. 
For the gravitino LSP, the produced axinos decay into gravitinos, yielding the hot dark 
matter. Thus, we obtain the constraint 

< 1.8 X 10-12 (i^) . (5.15) 

It is interesting to note that, if the axino is the LSP, the above cosmological constraints 
can be satisfied easily for Sq > lOHQeV. In addition, the axino can naturally explain 
the dark matter component of the Universe today. On the other hand, in the case of the 
gravitino LSP, the constraint is severer but not difficult to satisfy. For instance, models 
with Sq ~ lOiiQeV and ^3/2 ~ lOMeV will survive the constraint. The axino yield is 
further suppressed when Ba is less than 0.01. Here, one should note that the relation ( ^I^ ) 
leads to (ma/m3/2)5o < lO^^GeV for 5o < O.IX. 

Another process for the production of axino/gravitino is the LOSP decay. This becomes 
important either when Tu > Tj, or when < Tj and m„ > 2m-^. For 5*0 < O.IX, since x 
will dominantly decay into axinos with ^^^q ^ W~^T^^a, the gravitino production from 
its decay can safely be ignored. Let us examine the case that T^- is above Tf, for which 
LOSPs can be in thermal equilibrium with SM particles. The axino production is then 
determined by the total decay width of x 

r, . jj^. (5.16) 

A stau LOSP has r| = 0(10-^), whereas a bino LOSP has 0(10-1) < rf < 0(10-^), 
depending on their mass. The decay rate F^ is thus smaller than the Hubble parameter 
at r = To-, for So > lOi'^GeV. If the decay temperature is below m^, the yield of axinos is 
naively estimated as = + Y^, with Yq^-^ given by 



3 X 10-^ 



dx 



100 e-i/^ 



, (5.17) 



10-3 V200GeV/ \ 5n / 
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for X = T/rriy. with = Tfjjray^ and = Tjjra^. After freeze-out of ah the remaining 
LOSPs win decay into axinos. This contribution gives Yq ■, which is less than 0(10"^^) for 
Tj < m^/20. Hence, a constraint is placed on Y^. Substituting Y^ by Y^ in the equation 
(5.14), one can find that the LOSP decay would not cause cosmological problems for an 
axino LSP if < m^^/lO, or S^o > 10"GeV. For 5*0 ~ 10i°GeV, it is stih possible to 
achieve To- below m-^/10 when the saxion acquires mass, = 0(0. IMq). On the other 
hand, if the gravitino is the LSP, one would need < m;^/15 and 1713/2 <C IGeV to 
suppress the abundance of hot LSPs. Notice that the case with T^- > should obviously 
be excluded^^. 

Let's move on to the case that To- is lower than Ty, for which there are only a negligible 
number of LOSPs in thermal bath. However, LOSPs will be produced abundantly from the 
saxion decay if the decay process is kinematically allowed. In this case, the annihilation 
process can be effective to reduce their number density, depending on the decay width of 
the saxion into x- 

Po^^ = r^BaT^. (5.18) 

From (13), one obtains 0(10"^) < r| < 0(10"^) for a bino LOSP, while O{10-^) < r| < 
0(10^3) in the stau LOSP case. The LOSPs produced by the saxion decay will annihilate 
with each other if the interaction rate is much larger than F-j^. This condition translates 
into 

where (• • • ) represents the thermal average of the annihilation cross section times the 
relative velocity of x- For a bino or stau LOSP, the above condition is indeed satisfied 
well, implying that LOSPs are so abundant. Therefore, the annihilation process occurs 
quite effectively until the Hubble parameter becomes comparable to the annihilation rate. 
After annihilation, x decays into an axino and axion, and thus the axino abundance is 
determined by 

yx' ^Y^^\ ( 90 ^ 1 

^ 3 X 10-12 ^ . (5.20) 



For a stau LOSP, the annihilation cross section is roughly given by (cann'^reOf ~ lOa^m^^f ^ 
|3^, and has a value of 0{lQ-^)Ge\-'^ for nif = 200GeV. A bino LOSP has (crannWrei)^ ~ 



Aa{mimj'^ fS^, which is of C'(10-^)GeV-2 for mr = 500GeV. In addition, T„ cannot be 
much lower than O.lGeV for lO^^GeV < Sq < lO^^QgV. It is thus not difficult to make Y^' 
less than Y^ . This implies that, even for below Tj, the saxion can have mass ttt-o- > 2m^ 
without causing cosmological difficulties. 



Axinos can also be produced by thermal scattering |3^] after saxion decay. But the thermal production 
would be negligible when the saxion decay temperature is much lower than the MSSM sparticle masses. 
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We close this section by summarizing the implications of the Higgs coupling to S. As 
discussed already, the coupling between HuH^ and S naturally explains the fi and Bfi 
terms in the MSSM. Furthermore, it makes the saxion and axino interact with the MSSM 
particles via various couplings suppressed by the PQ scale. Mediated by these interactions, 
the LOSP decay into axinos can occur inside the detector, depending on the PQ scale. 
The decay length will give us a direct information on the PQ scale. On the other hand, 
the saxion products might cause cosmological problems, once the Universe is dominated 
by the saxion. Since the saxion decays into axions, axinos, and gravitinos, there are 
constraints from the Big Bang nucleosynthesis and the dark matter abundance. However, 
the cosmological constraints can naturally be satisfied when the jj, term is generated via 
the KN mechanism, for which the saxion coupling to the MSSM particles is stronger than 
those induced by the PQ messenger loops. 

6. Conclusion 

In this paper, we have studied a simple axion model that establishes a connection between 
the origin of the Higgs jJi/BjjL term and the solution to the strong CP problem within the 
framework of gauge mediation. Such a connection is possible if the model possesses the PQ 
symmetry under which the Higgs bilinear HuHd is charged. The PQ symmetry breaking 
is governed by SUSY breaking effects. We pointed out that a crucial role is played by the 
mixing between the messengers transmitting the SUSY breaking and the PQ symmetry 
breaking to the MSSM sector. In the presence of such mixing, the PQ scale is radiatively 
stabilized at a scale below the gauge mediation scale. This stabilization mechanism can 
apply to other CclS6S clS well, such as models with a generalized messenger sector, or flaton 
models where S corresponds to a flaton field. 

Also important is that the model provides a natural explanation for the presence of both 
H and Bjj,. They are generated with the correct size from a coupling between HuH^ and S, 
which also induces the saxion/axino interactions with the MSSM particles. Furthermore, 
the phase of B is aligned with that of the gaugino masses, thereby not spoiling the nice 
property of gauge mediation that the induced soft terms do not lead to excessive flavor 
and CP violations. In the model, the LSP is either the axino or the gravitino depending 
on the scale of gauge mediation, while the LOSP can be the bino or the stau. The Higgs 
coupling to S leads to that the LOSP mainly decays into axinos with coupling suppressed 
by the PQ scale. Thus, the collider signature highly depends on the PQ scale. If the saxion 
is stabilized at a scale around lO^^GeV or less, the LOSP can decay within the detector 
while giving distinct signals. On the other hand, for 5*0 larger than 10^*^GeV, the LOSP 
will decay with a rather long lifetime, but still a non-negligible amount of the LOSPs will 
decay inside the detector unless Sq is out of the axion window. 

We have also investigated the cosmological constraints placed on the saxion when 
it dominates the energy density of the Universe. In order for the saxion decay not to 
conflict with the successful predictions of the Big Bang nucleosynthesis, its branching ratio 
into axions should be suppressed. Moreover, the LSPs from the saxion decay should not 
overdose the Universe. In particular, the hot LSPs from the late decay of axino/gravitino 
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should be small enough to be consistent with the cosmological observations. All these 
constraints can be satisfied well for the PQ scale around lO^^GeV when the term is 
generated via the KN mechanism, i.e. from a superpotential Higgs coupling to S. This is 
because the saxion is coupled to the SM particles more strongly compared to the case of 
the GM mechanism. If there is an extra entropy production, models that incorporate the 
CM mechanism can still be cosmologically viable. 
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A. PQ sector soft terms 

In this appendix, we provide the expressions for PQ sector soft terms for the case that the 
saxion is stabilized at yiii\S\ < ©(A^/VStt^). The soft terms are parameterized in terms 
of {Mo, A^, N^, Aq,, Nq,}, and the supersymmetric couplings ga and y^^^. The Yukawa 
couplings are written as 

y^S^^ = ygSqq + yeSee, (A.l) 
for the PQ triplet q and doublet £. Integrating out $ + soft terms are generated as 



= iV, (lOiV^ + 1 + 16.^^) yl 
= {Id + - (5iV^ + l)yWs ) , 



"'^^'^^ ( lid - 2yl)9l + 4(5? - 2yl)9! + ^yl 



Ml " ' 15' 

= [lid - 2yM + ^{gl - 2yM +^yl], (A.2) 

at the scale just below A$, with 

e = {5mN^ + N^ + N^)yl + yl- (IGtt^)^!^, (A.3) 



\s\' 



where the gauge and Yukawa couplings are evaluated at A$ , neglecting the splitting in PQ 
Yukawa couplings. 
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At scales between A^, and A$, soft terms are RG evolved as 

"^^^ - ^ ^ {m^, + 2)ylA^ + 2N^yjA,-2(lglM^ + ^glMi 



dlnQ 87r2 {'"^'^ ' "'^t^^'^ ' -'-^^^^ - ' ^5 

dAf 



^ jsiVvpy^^, + {2N^ + 2)yjA, - 2 Qy^^s + ^fff^i) } 



dlnQ 

dlnQ 87r2 l^*? " V3 15^' ' 



dlnQ 87r2 [^'^^^ ^^3"^' ^' ' 10^ 

(3y2p- + 2y|P,-), (A.4) 



dlnQ 87r2 

for Pq = m| + m? + m| + | A^p, and = m| + m| + ml + \A£\'^. On the other hand, the 
running of PQ Yukawa couplings is determined by 



3N^y' + {2N^ + 2)yi - 2 -gi + -r^f , (A.5) 



dlnQ 87r2 ' ^^^'^ ' "^^^ ^ ^2^^ ' 10^ 

with gauge couplings given by 

1 ^ , /Avi>\ , ba , /MgutA 



- 2 + T-TT In ^ + — ^ In — — , (A.6) 



gliQ) 87r2 VQ/ 8^2 ^ Q 7 

at A^r < Q < A$. Here 6a are the beta function coefficients for the MSSM. 

B. Mixing parameters 

To derive the saxion/axino couplings to MSSM particles, one can make the replacements 
(3.13) and ( |3.14 ). The mixing between the saxion with neutral CP even Higgs bosons, h 



and H, is parameterized by N^u- 

{N^ , N^) = {—rifj sin a + n'^ cos q, ria cos a + sin a), (B-1) 
where rio- and are given by 

Ha = 2 sin(/3 — a) — j — r cos(/3 + a), 
1/^1 

= "'i' ~ (2 cos(/3 - a) - ^ sin(/3 + a)] , (B.2) 
mjj-mi V 1^1 / 

with a being the mixing angle for h and H. On the other hand, the parameters N?" for the 
axino mixing with the neutral gauginos and Higgsinos are determined by 

r.ra .rn x COS 2/3 [Mz . , Mz , \ 

(iV|o,7VL) - ^— ^(cos/?-nasin/3,sin/3-nacos/3), (B.3) 

"y. 1 — Ha Sm2p 
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for ria defined by 




Mz ( Mz 



Mz 



) 



(B.4) 



where Mz is the Z boson mass, and 6w is the weak mixing angle. Here we have neglected 
corrections suppressed by nia/Mj^yy or by maj ii-, and have used that there is mixing 
between the neutral Higgsinos and gauginos 



which arises after the electroweak symmetry breaking. 
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